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^ l' We derive the thermal correlators for twisted quantum fields on noncommutative spacetime. We 

, show that the thermal expectation value of the number operator is same as in commutative spacetime, 

but that higher correlators are sensitive to the noncommutativity parameters S^". 

00 , 

General arguments involving classical gravity and quantum uncertainties suggest that spacetime struc- 
ture should be "granular" at very short distances [T]. A specific model for this granularity is realized by 
the Groenewold-Moyal (GM) plane, where instead of the usual pointwise product (/ • g){x) on one 

Oh, works with the noncommutative product (/ * g){x) = /(a;)e^^'^^ ~^"g{x). A particularly important fea- 

ture of GM plane is that Poincare symmetries are automorphisms of the spacetime, albeit with a twisted 
coproduct A5i(A) [5] instead of the usual one. This in turn leads to deformation between the canonical 
(anti-)commutation relations [3l|4] in quantum field theory: 

>; Cpflq = rye^P^^aqflp, a^a^^ ■qe''P^'ia\^a[, 1 where p'' = (p°,p), p A q = p^Qf^^q^, 

■ ap^Jq = rje^^P^'^al^ap + {2n)^2po6'^{p — q) J and rj = ±1 for bosons/fermions. 

Op ' Since effects of noncommutativity become important at high energies, we expect that there may be im- 

portant implications in early cosmology with its attendant high temperatures. To this end, it is important 
CO ' to formulate the thermodynamics of such quantum field theories. Noncommutative spacetimes contribute 

. an additional subtlety to this issue, in that the usual facility of working with a finite volume V and then 

taking F — oo is not available to us. Thus the appropriate starting point for any discussion of quantum 
thermodynamics is the KMS condition (see for instance, [5Jin]). We shall demonstrate two different (but 
equivalent) methods of computing thermal correlators, based on dual versions of the KMS condition. 
Given an operator A (which may for instance be constructed from products of quantum fields, or from 
$H ' products of creation or annihilation operators) in the Heisenberg representation, its time evolution is given 

^ by A{t) = e^^'^Ae"'^'^, H = H — where H is the grand canonical Hamiltonian. It is important to 

emphasize that the r appearing in the above equation is not the coordinate time x'^, but the parameter of 
time evolution [7]- 

For any two operators A and B, we can define the retarded function 

GAB{r - r') EE ^^e{r - r')(([A(r), B{t')])) = -z0(r - r')[{A{r)B{T')) - 7y'(i?(r')A(r))] (2) 

where 9{x) is the Heavyside step- function. In situations where the Gibbs state p can be defined (for example 

for systems in a finite volume V), the thermal average of any operator X is {X) = ^'^'^ ^' = Tr [pX], 

where p = " ^" and Z — Tv [e~^^]. Advanced and causal functions can be defined similarly [5]. 

We will instead make use of the relation between Gab{t — t'), the thermal correlation functions 
^ab{t — t') = {A{t)B{t')) and Tba{t — t') = {B{t')A{t)), and the spectral density Jba{oj) defined by 
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JBA{^)e-"^^^~^'Uuj. (3) 
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Thermodynamic equilibrium (i.e. objects like Gab and Tab are functions of (r — r') only) and cyclicity 
of trace imply that 

/oo 
JBA{i^)e^'^e-'^^^-^'ULo, (4) 
-oo 

i.e. J-AB and J-ba satisfy the Fourier space version of the KMS condition 6 . We will use this as our 
starting point, rather than assume the existence of the Gibbs state p, thus circumventing the formal 
necessity of putting the system in a box of finite volume. 

For evaluating correlators of interest, we will follow the strategy outlined in [8]. The r- independent 
function 77' shall be chosen so that Gab satisfies a conveniently simple differential equation, as we shall 
show below. 

Heisenberg equations of motion for A{t) and B{t) imply that Gab satisfies 



.dGAB 



= 5{t r'){A{r)B{r) - r,' B{T)A{r)]) + (({A(r)H - HA{t)- B{r')})). (5) 



dT 

The Fourier transform Gab{E) = ^ G ab{t)c^^'^ dr can be written as 

Gab{E) ^ - / JsAiu^Ke^'^ "^'^IT^TTJ; 

using the integral representation 6{t — t') — ^ J^.^ — in ([4]). 

Using and the delta function representation S{x) — ^ | — STRe} S^* 

GABiu^ + ie) - GabIu^ - ie) - -zJi3A(c^)(e^'" - ??'), (7) 

which in turn gives 

TBAir-r') = r ^^""^^ ^ '1' ^^''^^ ' e-^(--')dig, (8) 

J-oo - V 

TABir'-T) = r Gab[E + »e) - Gab{E - ^e) ^-^E(r-r') ^j^^ (9) 
J-00 - V 

For a perfect quantum gas, the (grand canonical) Hamiltonian is 

(^^(^.-M)aU. (10) 

where the and at satisfy ([1]). Substituting j4(t) = Cp-^ (t), i?(T') = a]p^{T') in ([2]), we find that 

Gp.p. ^ -*^(r-r')[(ap,(r)at,,(r')) -??'(at.,(r')ap,(r))] (11) 

satisfies 

- {2nf2{p^^)5{T - t')5^{p^ - P2) + (a;p, - m) Gp,p,(r - r') (12) 

if we choose r]' — T^e"*^^^^^. 

The Fourier transform Gp-^^p^{E) of Gpjp2(T — r') is easily obtained: 



1 (2^)32(pio)(53(pi - P2 

= ■ 

Using ([5]) and putting r = r', we get 



G'pip2(£')-- r . (13) 

2tt E - (cjpj - p.) 



, t X (27r)32(pio)(53(p,-p3) 
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Since pi A p2 = if pi = P2, we have 



^ _ (27r)32(pio)j3(pi-p3) 



(15) 



which is same as the commutative correlation function. 

This result is not unexpected: translational invariance forces this upon us. Higher correlators however 
will not be so severely restricted by translational invariance. For instance, to calculate (ap^ap2'^P3^P4)j 
we substitute A{t) = ap^{T) and B{t') = ap^{T')a^p^{T')ap^{T') in ([2|): 

Gp.p,p.p3 = -^0(r-T') [(ap,(r)at^(r')at,^(r')ap3(r')) -??'(at,^(r')at^(r')ap3(r')ap,(r))] . (16) 
This satisfies 



^P4PlP2P3 



Sir - r')(27r)^ [2(pio)^^(p4 - Pi)(«L(T)ap3(r)) 



+ 27y(p2o)<5'(P4 -P2)e-*^^^^M«kW«P3W)] + K4 -M)Gp,p,p,p3 (17) 
for the choice r/' = rje-'P^^'-P^+P^-P^'l 



The Fourier transform Gp4pip2P3 is 



Gp.p,p.p3(£^) = ^;^^^^[2(pio)<53(Pi-P4)(<W«P3M^ 

+ 2,7(p2o)e^^^^^'-<53(p2-p4)(a^,(r)ap3(r))]. (18) 



Using ([U and putting t — t' we get 

/ t t \ (27r)3(2pio) (27r)3(2p2o) ^.3/ ^ 

flpi ap2 «P3 «P4 = — ^7 ^ 77 ^ [S Pi - P4 M'' P2 - P3 

+ 77e^P^^P^j3(pi-p3)^'(p2-P4)] . (19) 

This four-point correlator differs from its commutative counterpart by appearance of the ^-dependent 
phase e'^'i^P^ in the second term, and leads to interesting changes in observables like Hanbury-Brown and 
Twiss correlations 9 . 

Higher correlators may also be calculated by similar techniques as above. The computations are tedious 
but straightforward. Alternately, one can evaluate them by using the direct (as opposed to Fourier) 
formulation of the KMS condition |^. Let cj^.^t be a positive, linear, normalized map from the algebra of 
operators to C. For any two operators A and B we define two functions ^^^(t) and ^^b(''") 

GABir) = ^pAAr)B) - L.pAA)^pAB), (20) 

where, A{t) = e^^'^Ae"''^'^. The map w^^^ is a thermal state corresponding to the inverse temperature (3 
and chemical potential fi if 

+ */?). (21) 

Consider the operators a^. and a'^. , which stand for either creation or annihilation operators corresponding 
to the momentum state Pi and pj respectively. We define their twisted commutator as 

[a*,a*]e^a*a*~Ve^"-^-^^''^^P^a*a*, (22) 

where Q?p.p^ is 1 if a^. and a^. are of same type {i.e. are both creation or both annihilation operators), 
else is equal to —1. The commutation relations ([T]) imply that 



# # 
^Pi I '^pj 



|a^. ,a^. le — if "^pi and are of same type 



at. = Or 



= -77(27r)^2(pio)(5^(pi - Pj) if a# = a^^ and u,p. - u,p^. 
= (27r)32(pio)(53(p, - p,) if a* = ap, and a* - a^. (23) 
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Using (fTS]). we see that 



where 



\ait. , af 



' 1 - rye^-^Ki"^) 



a = 1 if a^. is a creation operator 

= —1 if a^. is an annihilation operator. (25) 

To compute the A^-point correlator ujj3_f^{a^^a^^ ■ ■ -flp^), we use ([22]) repeatedly to bring to the 
right side of the sequence. Linearity of i^p.^^ then gives us 

+ ^^-h^^^u;f,,,{alal--.a*^al) (26) 
where ^ on an operator denotes the absence of this operator from the sequence. The phase (jjj is given by 



^^ai.piApi. (27) 

For r = the KMS condition implies that 

ujpA^B) = Ljp^^{BAiil3)) = ijp^^{Be~^^Ae^'^). (28) 
For the Hamiltonian (|10l) . using (p4)) we can finally write 

N-l - 

^f^A«-0 = (E ■■•4+1 ■■•<j)C(/3,A^,^pJ, (29) 

where ^(/3, A^, Wp.) is given by 

C(/3,Af,Wp.) = ^ (30) 

This is the thermal version of Wick's theorem adapted to twisted quantum fields: the A^-point correlator 
is expressed in terms of the {N — 2)-point correlators. 
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